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ABSTRACT 

We present a detailed analysis of the intrinsic scatter in the integrated SZ effect - cluster mass 
{Y — M) relation, using semi-analytic and simulated cluster samples. Specifically, we investigate 
the impact on the Y — M relation of energy feedback, variations in the host halo concentration and 
substructure populations, and projection effects due to unresolved clusters along the line of sight (the 
SZ background). Furthermore, we investigate at what radius (or overdensity) one should measure the 
integrated SZE and define cluster mass so as to achieve the tightest possible scaling. We find that 
the measure of Y with the least scatter is always obtained within a smaller radius than that at which 
the mass is defined; e.g. for M200 (A^soo) the scatter is least for Fsoo (^1100)- The inclusion of energy 
feedback in the gas model significantly increases the intrinsic scatter in the Y — M relation due to 
larger variations in the gas mass fraction compared to models without feedback. We also find that 
variations in halo concentration for clusters of a given mass may partly explain why the integrated 
SZE provides a better mass proxy than the central decrement. Substructure is found to account for 
approximately 20% of the observed scatter in the Y — M relation. Above M200 = 2 x lQ^'^h~^yiQ, the 
SZ background does not significantly effect cluster mass measurements; below this mass, variations 
in the background signal reduce the optimal angular radius within which one should measure Y to 
achieve the tightest scaling with M2oo- 

Subject headings: cosmology: dark matter — galaxies: clusters: general — intergalactic medium — 
methods: N-body simulations 



1. INTRODUCTION 

The evolution of the number den sity of galaxy clusters 
is a sensitive co smological probe (iBahcall fc FanI 119981 : 
lEke et aT]ll998al ). As an indicator of the expansion rate 
as a function of time, the galaxy cluster number den- 
sity is sensitive to t he dark energy eq uation of state 
(|Haiman et all 120011: rVV^ller et al.ll2001h . This provides 
a growth-based dark energy test, an important comple- 
ment to the distance-based tests that have provided the 
most compelling eviden c e for dark ener gy to this point 
(jPerlmutter et all 119991: iSchmidt et all [1998). Galaxy 
clusters can be selected by many diverse methods, in- 
cluding (but not limited to) optical richness. X-ray ther- 
mal bremsstrahlung flux, and we ak lensing shear. The 
Sunyaev-Zeldovich (SZ ) effect (iSunvaev fc Zel'dovichl 
[197a iBirkinsha^ [19991 : iCarlstrom et all 120021) is being 
actively pursued as a method to detect galaxy clusters 
to high redshift, and e xperiments such as the Atacama 
Cosmology Telescope (Kosowsky 2003) and the South 
Pole Telescope (Ruhl 2004) are currently surveying the 
microwave sky to develop large catalogs of galaxy clus- 
ters that are uniformly selected. 

The key challenge for using galaxy clusters as precise 
cosmological probes is in understanding how to relate ob- 
servables (SZ flux. X-ray flux, optical richness, weak lens- 
ing shear, etc.) to a quantity that can be well predicted 
by theory, namely mass. The ultimate goal is to produce 
theoretical predictions of the distributions of observables 
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as a function of redshift and cosmological parameters. 
Short of this, one approach is to theoretically model 
the evolution of number density as a function of mass, 
and then estimate the mapping between observables and 
mass in order to predict the observed evolution. This 
mapping can either be estimated from theoretical consid- 
erations or be determined directly fr om the data, assum- 
ing some regularity in the mapp i ng ([Maiumdar fc Mohil 
l20mt lHul[2f)0l iLima fc Hij[2004l . [20051). In eithCT case, it 
is also important to understand the scatter in the m ap- 
ping between observable and mass ([Lima & Hu'2005). 

One reason that SZ surveys are an attractive possibility 
for dark energy studies is that the scatter in the mapping 
between SZ flux and mass is expected to be small from 
straightforward physical considerations ([Barbosa et al.l 
Il996l ): the integrated SZ flux is a direct indicator of the 
total thermal energy in the intracluster medium. It is 
exceedingly difficult to disrupt the global energy budget 
on this scale, as the total thermal energy of a galaxy clus- 
ter is on the order of lO^^^^'^ ergs. This expectation has 
been verified by hyd rodyna mical sim ulations of the SZ 
effect (SZE) ([Nagaj 2006; Motl et al.l 2005; White et al.1 
l2?)0l . 

It is important to identify sources of scatter; even 
though the scatter is expecte d to be small, it is still non- 
negligible l|Lima fc Hul 120051 ). Furthermore, the SZ flux 
integrated to infinity is not an observable, and the def- 
inition of the mass of a halo is ambiguous in a CDM 
paradigm (|Whitel[2001h . One must carefully define SZ 
flux and mass within some fiducial region. The mapping 
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between SZ flux and mass could plausibly be more or less 
tight (i.e. have more or less scatter) for different defini- 
tions of SZ flux and mass. In this paper we will address 
this question. 

The key physical processes that govern the density 
and temperature profiles of the hot baryonic gas in clus- 
ter interiors are still poorly understood and improving 
this s ituation i s currently a very act ive area of research 
(see iMotLsFal. 2 005: Bor gani et al. 2005; Nagai ,200^; 
Ettori et al.ll2006HS iiacki fc Springell [20061 : iRxmieo et all 
20061 iMuanwong et al., .2006[ and references therein) . 



Dark matter halos provide the potential well into which 
cluster baryons sink, setting the gas temperature scale. 
However, in order for theoretical work to be reconciled 
with the observations, it is clear that non-gravitational 
processes such as radiative cooling, star-formation and 
energy feedback (through supernovae, AGN outflows and 
galactic winds) play an important role (Evrard & Henryl 
1991; Kaiser 199 1;'Balogh et al. 1999; Eke e t al.. 1998B 
McCarthv et all 120031 : ISolanes et all I2005t iBode et'al.l 



20071) . Unfortunately, many of these processes involve 



small-scale, sub-grid physics that are very difficult to 
simulate directly. Typically a semi-analytic prescrip- 
tion must be incorporated in simulations to approximate 
these effects. 

Variations in the internal dynamics between clusters 
will be a significant contributor to the intrinsic scatter in 
the mass-SZ flux scaling relation. Furthermore, the large 
scatt e r in the halo mass -concentration r e lation (iJing 
2000; Eke ct al. 2001; Bull ock et all 120011 : iDolag et al 



2004 .Shaw et al...2006: Maccio et al.ll2007D and the wide 
range in substructure populations JPe Lucia et al.| [ ' 



iGao et al.l[200llGin et ail200l[Shaw'et al.ll2007n results 
in a gravitational framework that can vary signiflcantly 
between clusters of similar mass and redshift. In this 
study, we apply several different realisations of a semi- 
analytic model of intracluster gas, calibrated using obser- 
vations, to the output of a high-resolution N-body light- 
cone simulation to investigate some of these questions. 

Another significant contributor to the scatter in the 
SZ flux-mass scaling relation is errors in flux mea- 
surement caus ed by confusion due to pr ojection effects 
(jHolder et all 120071 ). iMotl et all (|2005D demonstrated 
that the thermal SZE integrated over a large fraction 
of the projected virial region of the cluster allows a ro- 
bust measurement of cluster mass, suppressing the im- 
pact of heating and cooling in the core. However, the 
noise due to background clusters will increase as one in- 
creases the angular size of the aperture within which the 
SZE is measured. It is important to assess the impact 
of these interlopers, particularly those clusters that are 
too dim to be identified individually and accounted for. 
In this study we investigate the angular scales at which 
scatter in the Y — M relation due to projection errors be- 
gins to dominate over that due to variations in internal 
cluster properties. 

The outline of this paper is as follows. In Section [2] we 
provide a brief overview of the mass-flux relation, and 
describe the simulations and gas model that we use to 
predict this relation. In Section [3] we investigate the in- 
trinsic scatter in this relation due to the choice of cluster 
mass definition, the area within which the SZE fiux is 
measured, the impact of variations in cluster gas temper- 
ature and density, and the variations in host halo con- 



centration and substructure populations. In Section [?] 
we use a simulated lightcone to investigate the impact 
of projection effects on mass estimates as a function of 
cluster and aperture angular size. Finally, in Section [5] 
we summarise and discuss our conclusions. 

2. METHOD 

2.1. Sunyaev-Zel'dovich Scaling Relations 

The thermal SZE is a distortion of the CMB caused by 
inverse Compton scattering of CMB photons (at temper- 
ature Tcmb) off electrons (at Tg) in the high temperature 
plasma within galaxy clusters. To first order, the tem- 
perature change at frequency of the CMB is given by 
AT^/T„„6 = Mx)y, where f^{x) = a;(coth(x/2) - 4), 
X — hv /kBTcmb, and y is the normal Compton parame- 
ter 
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n,{l)TS)dl , (1) 



where the integral is along the line of sight. The in- 
tegrated temperature distortion (the SZ flux) across the 
surface of a cluster is deflned as 
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where il is the solid angle subtended by the cluster on 
the sky. For a cluster at redshift z, 



Y. = 
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n,{l)T,{l)dV , 
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where dA{z) is the angular diameter distance, and the 
integral is now over the volume of the cluster. Yz is thus 
measured in angular units squared. However, it is often 
assumed that simulated clusters reside at a redshift of 
zero, and thus the l/(iyi(z)^ factor is omitted. In this 
case the units of Y are Mpc^ (henceforth represented by 
omitting the subscript z). In the following section, in 
which we investigate the intrinsic scatter in the Y ^ Al 
relation due to variations in internal cluster properties, 
we adopt the latter definition. In the final section, where 
we investigate the impact of projection effects on the 
Yz — M relation using a simulated lightcone, we measure 
the integrated Compton parameter as defined in Eqn. [3l 
The mass of a halo is typically defined as the mass con- 
tained within a region of spherical overdensity A times 
greater than the critical density. 

Ma = ^7ri?iApc(z) , (4) 

where Pc{z) is the critical density at redshift z. 

From spherical collapse theory, the region in which 
matter is predicted to have virialized is defined at z = 
by the overdensity contour Ac = 178 in a D,m — 1 



universe or Ac 



96 in a flat 
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0.26 universe 



(iLahav et al.l fl99H iLacev fc Cold Il99l iCole fc LacevI 



ll996l:lEke et al.lll996l : iBrvan fc Normanlll99 8^. Most re- 
cent studies of halos in simulations deflne mass adopting 
either the redshift dependent Aei or a flxed overdensity 
factor of A = 200. lEvrard et al] (|2007D compare common 
definitions of halo mass using a suite of N-body simula- 
tions, focusing specifically on the the M — ctdm virial 
relation, where ctdm is the one-dimension dark matter 
velocity dispersion. They constrain the region of mini- 
mum variance in the virial relation to being within the 
overdensity contour A = 200. 
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In the absence of non-thermal heating and cooling pro- 
cesses, we expect a self-similar scaling between cluster 
mass and SZ flux: 

Yo,f,,M'/'^E{zr/\ (5) 
where /gas is the cluster gas fraction, and E{z) = 
(r2M(l + z)"^ + J^a)^^^- In practice, many studies of simu- 
lated clusters have measured a steeper slope than (5/3), 
attributed to the presence of radiative cooling and non- 
gravitational heating in t heir simulations (IWhite et ahl 
l200llda Silva et al.li2004l: iMotl et al.l[2005l : lNagaill2006ir 

To date, there have been only a limited number of clus- 
ters o bserved using the SZE. Recently, iMorandi et al.l 
in an X-ray and SZ study of 24 clusters in the red- 
shift range 0.14-0.82, found a good agreement between 
the slope of the observed Y — M relation and the self- 
similar predictions for their cooling core sample. How- 
ever, they find a significantly shallower slope when they 
include both cooling-core a nd non cooling-core clusters 
in their analysis (see also iBenson et al.l 1200^ . They 
also measure greater scatter in the Y — M relation than 
is typically seen in simulations. 

The existence of a tight SZ flux-mass scaling relation 
provides the crucial link between SZ cluster surveys and 
constraining cosmological parameters. It is therefore im- 
portant to understand the character and physical origins 
of the scatter in this relationship, and how one can sup- 
press its impact to obtain an accurate measure of cluster 
masses. In particular, there has been little investigation 
concerning at what radius (or equivalently, what over- 
density) one should attempt to measure the integrated 
SZ flux, and define cluster mass, so as to achieve the 
tightest possible scaling. We investigate this issue in de- 
tail in Section [31 

2.2. Simulations and Gas Model 

For this study we adopt two complementary ap- 
proaches to measuring and characterising the scatter in 
the Ya — scaling relation. Our main cluster sam- 
ples are generated through applying different rcalisa- 
tions of an intra cluster gas model ( Ostrikcr et al. 2005, : 
iBode et al]|2007[ ) to individual dark matter halos iden- 
tified in a high resolution N-body simulation. For com- 
parison, we also use a second sample of clusters taken 
from the output of an SPH simulation. In this section 
we describe in more detail the simulations and gas model 
used in this analysis, and the derived cluster catalogs. 

To create our main simulated cluster catalog, we begin 
with the output of a large {N = 1024^ particles) cos- 
mological dark matter simulation. The cosmology was 
chosen to be consistent with that measured from the 
3rd- year WMAP d ata combined with large-scale struc- 
ture observations ([Spergel et al.l [20071 ) . namely a spa- 
tially flat LCDM model with parameters: baryon density 
J7b = 0.044; total matter density flm = 0.26; cosmolog- 
ical constant JIa — 0.74; linear matter power spectrum 
amplitude erg = 0.77; primordial scalar spectral index 
Us = 0.95; and Hubble constant Hg = 72kms~^Mpc~^ 
(i.e. h = 0.72 = Ho/lOOkms^^Mpc-i). The simulated 
volume is a periodic cube of size L = 320ft.~^Mpc; the 
particle mass nip = 2.2 x 10^h~^MQ, and the cubic spline 
softening length e = 3.2h~^kpc. 

The matter distribution in a light cone covering one 
octant of the sky extending to z — 0.5 was saved in 



315 time slices. Dark matter halos were identified using 
the Friends-of-Friends algorithm with a comoving link- 
ing length parameter b — 0.2. For analysis presented 
here, we select a sample containing the 1267 lowest red- 
shift clusters of mass M^r > 5 x 10^^h~^MQ, spanning 
a redshift range of 0.002 < z < 0.14. 

The cluster gas distribution in each halo was calculated 
using the semi-analytic model described in lOstriker et al.l 
(|2005) and Bode et al. (2007). In brief, a 3D mesh 
(with cell side-length 4e = 12.8h~^ kpc) is placed around 
each cluster, with the gas pressure and density deter- 
mined in each mesh cell assuming hydrostatic equilib- 
rium and a polytropic equation of state (with adiabatic 
index F = 1.2, [Ascasibar et al. 2006). It is also assumed 
that the gas in the densest cluster regions has cooled and 
condensed, formi ng stars. At z = 0, the stellar/gas mass 
ratio is set to 0.1 ()Lin et al.ll200llVoevodkin fc Vikhlillh] 
I2004D . To compute the star/gas ratio at z > 0, the star 
formation rate was as sumed to follow a delay ed exponen- 
tial model (Eqn. 1 of lNagamine et alll2006l ). with decay 
time T = 1.5Gyr. 

As discussed in detail in iBode et al.l (|2007t ) , the most 
important free parameter in this model is the energy in- 
put into the cluster gas via nonthermal feedback pro- 
cesses, such as AGN outflows and supernovae. This is 
set through the parameter ef, such that the feedback 
energy is efM^c^, where M* is the stellar mass in the 
cluster. In this study, we set e/ = 4 x 10~^, which is a 
value that provides good agreement between simulation 
cluster prope rties and observe d X-ray scaling relations 
(see Fig. 2 of IBode et aL[[2007[) . 

To investigate the impact of the feedback parameter 
and the assumption of hydrostatic equilibrium on the 
Y — M scaling relation, we have created three realisations 
of our cluster sample. In the first, the gas distribution in 
each halo is calculated using the full model (for which the 
measured X-ray scaling relations agree well with obser- 
vations). For the second sample we remove the effects of 
feedback by setting £/ to zero, so that the gas is merely 
re-arranged into hydrostatic equilibrium having removed 
the fraction that is assumed to have formed stars. 

For our final sample, we take an very simplistic ap- 
proach, assuming an isothermal model in which the gas 
directly follows the dark matter distribution. This is 
achieved through setting the gas density in each mesh 
cell to rif, / riuAf times the interpolated dark matter den- 
sity. The gas temperature is constant with radius, with 
the global temperature set according to the scaling rela- 
tion measured by Vikhlinin et al. (2006), 



E{z)M5aQh-^yiQ = A 



X 



kT \ ' 

f^-^ raw \ 
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where Ax — 3.32 x lO^'^ft-^^Mo, T,nw is the gas mass 
weighted temperature, and the slope a — 1A7. We em- 
phasize that this final sample does not represent a phys- 
ically accurate model, and is included as an aid to un- 
derstand the results of the more sophisticated approach. 
Hereafter, we refer to each sample as FULL, NOFB and 
BASIC, respectively. 

In practice, Y is measured for the clusters in each of 
our samples by summing up the electron pressure UeTe 
in each mesh cell along one direction, thus creating a 
2-d y image of the cluster (for example, see Fig. 7 of 
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lOstriker et al.1 [20051) . Note that this is directly equiv- 
alent to an integral along the line-of-sight through the 
cluster. We then calculate the integrated Y parameter, 
by summing all cells that lie within circular aperture of 
radius Ra (defined below) from the cluster centre. 

2.3. Hydrosimulation 

For comparison with the gas model results, we also 
make use of a (publicly available^) sample of clusters 
from the output of an smoothed-particle-hydrodynamics 
(SPH) simulatio n , evo lv ed using the HYDRA cod e 
ICouchman et all (|1995f ): iPearce fc CouchmarJ (|1997D . 
The simulation contained 160'^ gas and 160^ dark mat- 
ter particles in a box of side length 100/i-iMpc. The 
cosmological parameters were J^m = 0.35, fif, = 0.038, 
flA = 0.65, (78 = 0.9 and h = 0.71. The gravitational 
softening length was set to 25h~^ kpc, once the simula- 
tion had evolved past a re dshift of 1. The simulat ion is 
described in more detail in iMuanwong et all (|2001h . 

It is important to note that this simulation did not 
contain any radiative cooling and so the dense, cold gas, 
in the centre of clusters is unable to cond ense and form 
stars,^ IMuanwong et all (|2001l . |2002[ ) and iThomas et all 
()2002| ) find that the inclusion of radiative cooling removes 
this gas and heats the surrounding material, providing a 
temperature and entropy boost to the ICM. Inclusion of 
radiative cooling was also found to be necessary in order 
to match observed X-ray scaling relations. Therefore, the 
cluster gas in the SPH simulation used here is likely to 
be colder and significantly more centrally concentrated 
than expected in more physically realistic models. 

Cluster particles are selected by identifying all par- 
ticles within a cube of side-length ARyir around each 
cluster centre in the z = simulation snapshot, where 
the virial radius a nd cluster centre a re tak en from the 
cluster catalog of IMuanwong et al.) ()200 1'). In total, 
our sample contains 212 clusters with a minimum virial 
mass of 2.5 x 10^^h~^MQ (where the virial overdensity 
« 108). All clusters contain at least 1000 particles 
(dark matter plus gas) within their virial radius. 

To measure the integrated SZ flux we place a 3d mesh 
around each cluster, calculating the interpolated gas den- 
sity and temperature wit hin each cell according to th e 
smoothing kernel given in IThomas fc CouchmanI (|1992f ). 
For consistency we set the mesh cell side-length equal to 
that used in the gas model described in the previous sec- 
tion (12.8/i~^ kpc). The integrated SZ Y parameter is 
then calculated for each cluster from the mesh in exactly 
the same way as for the gas model samples. Henceforth, 
we label this sample SPH. 

3. RESULTS 

3.1. Measuring Y — M scaling relations 

For each cluster sample, we measure both the mass Ma 
and the integrated Compton Y parameter Ya within re- 
gions of density Apc, where the overdensity varies within 
the range 50 < A < 1500 - corresponding to roughly 
0.37 < r/i?200 < 1-77 - in steps of A = 100. Thus for 
each cluster we obtain 16 measures each of M and Y. 
From these various measures we will determine empiri- 
cally the combination of Ma and Ya that provides the 

^ taken from |http : //vlrgo . susx. ac .uk/clusdata. html | 



lowest scatter around a power-law relation. The outer- 
most value of A = 50 is outside the virial overdensity 
calculated from spherical collapse, Ac « 95. The in- 
nermost value of A = 1500 was chosen for two reasons. 
First, we are constrained by the size of the mesh cells in 
the gas model. Second, none of our cluster samples are 
expected to correctly reproduce the gas temperature and 
density in the cluster core (r < 0.1i?2oo), due to the com- 
plicated interplay between cooling and heating processes 
in these regions. Hence we pick an inner boundary which 
is considerably outside the cluster core. 

For each of the 256 combinations of Ay and Am, we 
fit their relation with a power-law. 



Y^E{zf^HO^ 



Ma 
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where 10'^ is the normalisation at lO^'*ft-~^M0 and a is 
the slope. We define the scatter around the best fit rela- 
tion as; 
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where Yay is the fiux measured within Ray^ ^Am is 
the fitted fiux for cluster of mass Mam N is the 
total number of clusters in the sample. The fitting is 
performed using a Levenberg-Marquardt algorithm. To 
avoid introducing a bias into our results by having de- 
fined the minimum cluster mass in our sample at a partic- 
ular overdensity, we fit (and calculate uym) using only 
the 1000 most massive clusters at each Am in our gas 
code samples. This corresponds to a limiting mass of 
5 X IO^^/i-^Mq at Am = 200. For the SPH sample, we 
take the 150 most massive clusters, corresponding to a 
minimum mass of 3 x IO^^H'^Mq at Am = 200. 

We note that for all combinations of Ay and Am, a 
power-law is found to be a very good fit to the Y ~ Al 
relation. The scatter, S = \nY — InY , around this re- 
lation in the vertical direction is gaussian with a slight 
tail towards low (high) values of Y for low (high) Am- 
Hence the scatter, CTyM, in Eqn. [8] represents the stan- 
dard deviation of the distribution of IuYay around the 
power-law fit. 

3.2. Slope and Normalisation 

To enable comparison of the results of our gas model 
and SPH simulation clusters with previous work, in Table 
[T] we give the slope a and normalisation index A (and 
their associated errors) from the power- law fit to the Y — 
M relation at selected values of Ay = Am- 

At A = 200, the slope for the FULL cluster sample 
(1.70) is close to the 5/3 expected from spherical collapse 
theory. Within smaller radii, the slope becomes steeper. 
This is due to the increasing impact of energy feedback 
on the gas fraction, /gas = Mg(< R)/Mtot{< R) within 
decreasing R. The additional energy increases the gas 
temperature, inflating the overall gas distribution and 
decreasing the baryon fraction within any given radius. 
For lower mass clusters this effect is more pronounced, 
thus ste epening the slope of the Y — M relation (see Fig- 
ure 2 of lBode et al.ll2007h . The slope a becomes greater 
than 5/3 due to the weak mass dependence of the gas 
mass fraction. This mass dependence has been observed 
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50 


1.51 (±0.05) 


-5.59 (±0.01) 


1.49 (±0.03) 


-5.52 (±0.01) 


1.64 (±0.02) 


-5.59 (±0.01) 


1.64 (±0.05) 


-5.71 (±0.01) 


100 


1.62 (±0.05) 


-5.64 (±0.01) 


1.51 (±0.03) 


-5.50 (±0.01) 


1.64 (±0.02) 


-5.54 (±0.01) 


1.64 (±0.05) 


-5.63 (±0.01) 


200 


1.70 (±0.05) 


-5.62 (±0.01) 


1.51 (±0.04) 


-5.44 (±0.01) 


1.65 (±0.02) 


-5.46 (±0.01) 


1.67 (±0.06) 


-5.54 (±0.02) 


500 


1.83 (±0.05) 


-5.50 (±0.01) 


1.53 (±0.03) 


-5.30 (±0.01) 


1.66 (±0.02) 


-5.29 (±0.01) 


1.72 (±0.05) 


-5.39 (±0.02) 


1000 


1.92 (±0.04) 


-5.36 (±0.01) 


1.55 (±0.03) 


-5.19 (±0.01) 


1.65 (±0.02) 


-5.16 (±0.01) 


1.77 (±0.04) 


-5.25 (±0.02) 


1500 


1.98 (±0.04) 


-5.26 (±0.02) 


1.57 (±0.03) 


-5.12 (±0.01) 


1.56 (±0.03) 


-5.10 (±0.01) 


1.77 (±0.04) 


-5.17 (±0.03) 



TABLE 1 

Slope a and normalisation index A measured for a selection of overdensities for each cluster sample. 



in hydrodynamical si r nulations (|A /[uanwo n g ; et al.l (|2002l ) ; 
da Silva et all (l2004D :l Kravtsoy et al. (20051) t hough see 
Ettori et al.l (|2006[ )'). We refer the reader to iBode et all 
( 2007t l for a more extensive discussion on the gas frac- 
tion obtained from this model, and comparisons to other 
studies and observations. 

In general, both the slope and normalisation that we 
measure for this m odel ar e in rea sonably good agreement 
with the results of iNagail ()2006[ ) for their hydrodynamic 
simulations incl uding cooling, star fo r matio n, and super- 
nova feedback (jKravtsov et alJ 12001 120051 ) . differing by 
less than 6% in a and less than 1% in A at A = 200 and 
500. Ida Silva et al.l (|2004D investigated the impact of in- 
cluding radiative cooling or pre-heating the gas in their 
simulations on the 1200 ~ -^^200 relation, finding slopes 
of 1.69, 1.79 and 1.93 for their adiabatic, cooling and 
preheating runs respectively (their adiabatic simulation 
is the same as that used in this study). They also find 
that the increased slope is due to a decrease in the gas 
fraction in lower mass (M200 < 10^'^h~^M(7)) clusters in 
the the cooling and preheating simulations. The steeper 
slopes measured compared to this study are probably due 
in part to the lo wer mass th r eshold 10^^h~^MQ in their 
cluster samples. iMotl et aD ()2005l ) measure an increase 
of 0.1 in slope between their adiabatic and cooling sim- 
ulations, although the slope decreased again once they 
allowed for star- formation in their simulations. 

The slopes measured from the NOFB sample are typ- 
ically much lower than those of our other samples or in 
previous studies, indicating that assumptions of hydro- 
static equilibrium alone (i.e. without feedback) do not 
reproduce the results of more sophisticated simulations. 
Indeed, more realistic slopes are measured from our BA- 
SIC sample, although this is likely due to the way in 
which cluster temperature in this sample is calibrated 
from observational results. 

The results for our SPH sample match closely those 
of previous adiabatic simulations, although we note the 
lower normalisation measured her e compared to the adi- 
abatic simulations of lNagail (|2006l ) {A = —5.54 compared 
to -5.42 at -R200: corresponding to a decrease of 25% in 
Y); the difference is probably due to the different ratios 
of rib/ilM between the simulations (0.11 and 0.14, re- 
spectively). This demonstrates the sensitive dependence 
of the normalisation of the Y — M relation on fib/f^M- 



3.3. Scatter inY — M relation 

Figure [T] displays the scatter ctym hi the Y ~ M re- 
lation obtained for all combinations of Ay and Am, for 
each of our cluster samples; FULL {top left), NOFB {top 
right), BASIC {bottom left) and SPH {bottom right). We 
consider first the results of the BASIC sample as it cor- 
responds to the most simplistic case. As described in 
Sec I2.2[ for this sample we make the simple assumption 
that the gas in clusters is isothermal and follows the dark 
matter density. Therefore, variations in the dark matter 
mass within Am are also reflected in the measured flux 
within the same overdensity, explaining the valley in the 
Am- Ay plane along Am = Ay- This characteristic is ev- 
ident in the results for all four of our cluster samples. The 
minimum scatter for the BASIC sample is near A = 500 
for both mass and flux. This is because the tempera- 
ture for these clusters was set via the M — Tmw scaling 
relation, with the mass defined at A = 500. 

The NOFB sample exhibits similar features to the BA- 
SIC sample, the most significant difference being the 
smaller rise in scatter as Ay increases at low values 
of Am- This is straightforward to understand: assum- 
ing hydrostatic equilibrium results in gas density profiles 
that are less sensitive to variations in the internal mass 
distribution of the host dark matter halo, especially in 
the inner regions, than for the clusters in the BASIC 
sample. The gas density and temperature profiles are 
set by a combination of the total gravitational and ki- 
netic energy of the gas and the surface pressure at the 
virial radius. We find that our recipe for star- formation 
has little effect on ctym- 

The FULL cluster sample - for which we added a pre- 
scription for energy feedback in the gas model - resembles 
the NOFB sample, but with uniformly greater scatter. 
Including feedback results in more extended (or 'puffed 
out') gas distributions, \ yith a shallower g as density pro- 
file in the inner regions (jBode et al.ll2007l ). Importantly, 
we have found that this substantially increased scatter is 
due to the greater variance in the gas fraction within any 
given overdensity. Table [2] gives the mean (/gas) and frac- 
tional standard deviation (7gas/(/gas) of the gas fraction 
within selected overdensities for the FULL, NOFB and 
SPH samples. Due to the dependence of /gas on cluster 
mass in the former, we calculate these quantities for clus- 
ters in the mass range lO" < Af2oo < 2 x lO^'^/i'^Mg. 
Note that the mean values we measure within i?5oo and 
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Fig. 1. — The scatter, ctym (Eqn. [Sjl, around the mass - integrated SZE relation for cluster mass measured within the radius defined by 
the overdensity factor Am and the integrated SZE measured within projected radius defined by Ay. Each panel gives the results for one 
of our cluster samples: FULL (top left), NOFB (top right), BASIC (lower left) and SPH (lower right). 

mass is responsible for the increase in ctym over the three 
models in Table [21 

In the bottom right panel of Figure [T] we plot the (Tym 
plane for the clusters extracted from the SPH simula- 
tion. The results are very similar to those of the BASIC 
cluster sample. The scatter is least along the valley de- 
fined by Am = Ay, rising to a maximu m oi a ^ 0.4 on 
either side. As discussed in Sec. [2l Pear ce et al.l (|2000), 
and iMuanwong et al.l (|2001l . 12003 ). the omission of ra- 
diative cooling in this simulation results a large central 
concentration of baryonic gas (compared to models that 
allow for cooling and star-formation) and thus in a ICM 
radial density profile t hat closely resembles that of the 
dark matter (see also, Ida Silva et al.|[200l : lEttori et al.1 
|2006[) . Both the BASIC and SPH samples consist of clus- 
ters with very high gas density in their centres; this is 
not the case for the clusters in the other two samples. 

In the upper panel of Figure [2] we plot the overdensity, 
Aymin, at which the least-scatter measure of the inte- 
grated SZ flux is obtained for a given Am- In the lower 
panel we plot the corresponding value of ctym • The solid, 
dotted, dashed and dot-dashed lines represent the FULL, 
NOFB, BASIC and SPH cluster samples, respectively. 
For all four samples Avmin increases with Am, with a 
slope greater than one. For the FULL sample, the best 
measures of M50, Af200, and M500 are ^300, >5oo, and 
^1100) respectively. For Am > 700, Aymin appears to in- 
crease beyond the scales probed here. Performing a linear 
fit to this point gives Aymin = 1.9Am + 150. Overall, our 



TABLE 2 

Mean and fractional standard deviations of the cluster 
gas fraction measured within four radii. 



.R ino(~ Rv\v) are i n good agreement with the simulations 
of iKravtsov et al.l (|2005D . ahhough we measure greater 
scatter. For all three samples shown, the fractional scat- 
ter in /gas increases with overdensity, most noticeably 
for the FULL sample. At all radii, the percentage scat- 
ter in /gas for the FULL clusters is approximately 1.5 
times greater than that for the NOFB sample and more 
than 3 times than that of the SPH sample. We find that 
Cgas/ (/gas) increases when recalculated using lower mass 
clusters. Repeating this exercise for gas temperature, we 
find the scatter to be similar for all three models. Hence, 
we conclude that the scatter in /gas for constant cluster 
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Fig. 2. — (upper) The overdensity AYmin for which the scatter 
(Tym is least for a given Am- (lower) the corresponding value of 
(JYM for each combination of Am, AYmin- Solid, dotted, dashed 
and dot-dashed lines correspond to the FULL, NOFB, BASIC and 
SPH samples respectively. The light-dashed lines are the results 
obtained when only clusters within a narro w co ncentration range 
are selected from the FULL sample (see Sec. 13.41 1. Note that where 
^Ymin appears to become greater than 1500, the value of (tym in 
the lower panel is calculated at this overdensity. 

results show that the best measure of Ma is obtained by 
measuring Y within a smaller area than i?A • 

The lower panel indicates that for the FULL and 
NOFB models, the scatter is least for Am = 50, and 
is roughly constant for Am > 200. As discussed above, 
much of the scatter in these models is due to variations 
in the baryon fraction within each radius. This source of 
intrinsic scatter is suppressed in the very outer regions 
of clusters as the baryon fraction approaches the cosmic 
mean. For the SPH sample, the scatter decreases with 
increasing Am to a constant value of 15%. For the BA- 
SIC sample, the dip in ctym at around Am,y = 400 (as 
noted above) is clearly evident. 

3.4. Impact of halo concentration 

The concentration of a dark matter halo is normally 
defined as the ratio of the virial radius to the NFW scale 
radius, Rvir/fs and is thus a measure of the d ensity in 
the halo core regions avarro et al.lll996l[l997l) . In their 
early studies, NFW postulated that the concentration of 
a halo is an indicator of the me an density of the universe 
at the time of its collapse (Nav arro et al.l[r997[ ). Hence, 
present-day clusters with a high concentration will have 
formed early and remained undisturbed by major merg- 
ers, growing through gradual accretion and accumulation 
of much smaller objects. Halos with a low concentration 
will have had a more tumultuous recent merger history, 
and may still be in the process of relaxi ng i nto dynam- 
ical e quilibrium (IWechsler et al.ll2002, . 20061 : IShaw et all 
[200llMacci6 et al.ll2007D ! 

There have now been many studies of distribution 
of hal o concentration and its correla t ion with cluster 
mass (iJing 2000; Bull ock et aP 2001; Eke et all 120011 : 
iDolag et al. 2004; Avila-Reese et a l. 20051). Most stud- 
ies find that at a given redshift the distribution of halo 
concentrations is well described by a log-normal func- 
tion of dispersion w 0.22 and mean value < c 5 at 



z=0 (Ijh^l2000l : iBullock et al.l[200l . Halo concentra- 
tion is found to decrease as a function of increasing halo 
mass and is well describ ed by a power-law of slope in the 
range -0.14 to -0. 1 (Bullock et al.ll20()it iDolag et al.ll2004l : 
IShaw et al.l[200llMacci6 et al.ll2007D . However, there is 
typically much scatter in this relation. One might ex- 
pect to find that variations in cluster concentration (and 
thus central potential) will have a significant effect on 
the gas temperature and density profiles, and thus the 
Y — M scaling relation. We therefore fit NFW profiles 
and measure the concentrations of all the halos in our 
samples, in order to quantitatively measure the impact 
of this variation. 

In Fig. 3, we recalculate ctym for the FULL cluster 
sample, using only clusters with concentration 5.25 < 
c < 5.75. The clusters in this sample are thus self-similar 
in terms of their dark matter mass distribution. By com- 
paring this figure with the upper-left panel of Fig. [H it 
is clear that constraining halo concentration significantly 
changes the geometry of the ctym plane. For constant 
Am, the scatter decreases with increasing Ay- For con- 
stant Ay, cTyM increases slowly as Am decreases towards 
Am = 500, and then falls off rapidly in the outer regions. 
We have verified that similar results are obtained for dif- 
ferent concentration bins and our three other cluster sam- 
ples. On removing clusters from our sample for which the 
NFW profile is a poor fit, the peak at Am = 500 disap- 
pears, hence this feature seems to be related to clusters 
for which we have a poor measure of concentration. 

The light-dashed line in the upper and lower panels 
of Figure [2] give Ay,„in and the correspond scatter at 
each Am for the clusters in this concentration range. 
For all but the two outermost radii, the optimal radius 
within which to measure the integrated SZE appears to 
be less than i?i5oo- (Note that where Aymin appears 
to become greater than 1500, the value of gyM in the 
lower panel is calculated at this overdensity). iMotl et ahl 
(|2005f ) demonstrated empirically that the integrated SZE 
(within i?5oo) provides a more accurate measure of clus- 
ter mass than the central decrement, y. Our results sug- 
gest that variations in halo concentration (and thus cen- 
tral potential) between clusters of the same mass may 
be respon sible for mu c h of t he very large scattered ob- 
served by IMotl et al.] ()2005f ) in the y-M relation. In- 
deed, they find that I500 is less sensitive to mergers than 
y; it has been previously shown that concentration is 
strongly influenced by the dynamical state of a halo (e.g. 
IShaw et al.ll200fil : iMaccio et al.ll200l . 

3.5. Impact of Substructure 

Halo substructure i s iden tified using the SKID algo- 
rithm of lStadel et al.l ()1997l ). with a smoothing length of 
5e, where e is the spline kernel force softening length of 
our simulation. The minimum mass of subhalos that we 
resolve are « 3 x 1O^^M0, corresponding to 100 simula- 
tion particles. For each halo, we define the substructure 
fraction fs = Mg/Myir: where Ms is the total subhalo 
mass. We find a mean value fs = 0.084 for our N-body 
halo sample, where we have included all subhalos with 
centres located within the radius i?5o of the cluster. This 
is in good agreement with previous stu dies of the sub- 
structure content on cluster mass halos (iDe Lucia et alJ 
l200llrTao et al.ll200llC^ll et al.ll2004l : IShaw et al.ll200l 7 

We find that, on removing all clusters with fs > 0.084 
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Fig. 3. — The ctym plane calculated for the FULL cluster sample 
having included only clusters with halo concentration in the range 
5.25 < c < 5.75. It is clear that the geometry of the plane is 
significantly altered compared to the results when all clusters are 
included, as shown in (the upper-left panel of) Fig. [T] 

from the FULL sample, (Tym is uniformly reduced by 
10%. Taking the sample quartile with the lowest sub- 
structure fraction (/s < 0.050) we find the average de- 
crease in (TYM is 20% (e.g. at Ay = Am = 100, (tym 
is reduced from 0.30 to 0.25). Taking the quartile with 
the highest substructure fraction (/^ > 0.010) increases 
cym by 32% (compared to the result for the entire sam- 
ple). Note that in all cases the geomentry of the <jym 
plane does not change significantly from that shown in 
the upper-left panel of Fig. [1] 

Overall, wc find that substructure accounts for approx- 
imately a fifth of the total amount of scatter in the Y — M 
relation, but - unlike the variations in halo concentration 
- does not strongly effect the dependence of ctym on the 
combination of Ay and Am. 

4. IMPACT OF PROJECTION EFFECTS ON THE Y - M 
RELATION 

Taking the standard definition of cluster mass, M200, 
we have established that the intrinsic scatter in the 
Y — M200 scaling relation is least when Y is measured 
within i?5oo (~ 0.66R2oo)- Assuming that ctym does not 
change significantly with redshift, for a cluster at red- 
shift z this radius translates to an optimal angular ra- 
dius, Oopt^ within which one can obtain the least-scatter 
estimate of cluster mass through an SZE measurement. 
However in order to make a useful estimate of 9opt, we 
must also account for the impact of the SZ background - 
the superposition of faint fore- and background clusters 
along the line of sight - when measuring the integrated 
SZ flux of a cluster withi n some angular radius 9. 

Recently, iHolder et al.l ()2007f ) investigated the confu- 
sion in cluster SZE flux measurements due to the SZ 
background using a sample of sky maps constructed to 
represent a range in as of 0.6-1. Intr acluster gas was 
assumed to follow the analytic model of 'McCarth v et aTl 
(|2003 ). They found that the mass scale below which the 
rms fractional errors in flux measurements become less 
than 20% occurs at just above IO^'^H'^Mq, although this 
is sensitively dependent on as and increases with decreas- 



ing redshift. iHallman et al.l (|2007D explored the contri- 
bution of both gas in low mass (< 5 x lO^'^M©) halos and 
gas outside of cluster environments on the SZE signature 
of resolvable clusters, using lightcones constructed from 
an adiabatic hydro simulation. They find that the inte- 
grated background SZE makes up between 4% and 12% 
of the total cluster signal, depending on the beam size 
and sensitivity of a survey, although these values are av- 
eraged over a range of cluster masses. Both these studies 
concentrated on the SZ background as the main source of 
scatter in SZ cluster mass measurements. In this Section, 
we combine the intrinsic scatter in the Y — M relation 
due to internal variations between clusters (ctym) with 
that due to confusion with the SZ background, in order 
to obtain the optimal angular radius Oopt within which 
the integrated SZE can best be obtained for a cluster of 
given mass and redshift. Henceforth, we fix our definition 
of cluster mass to the standard definition, M2oo- 

For this purpose, we use the redshift range < z < 1.5 
of a lightcone constructed from an N-body simulation 
with the same cosmology as that used in the previous 
section, but with lower mass resolution (particle mass 
rup = 1.22 X 10^^ /i~^M0). The lower resolution was 
enforced due to the amount of disk space required to store 
the lightcone over an octant (w 5250 square degrees) 
out to high redshift, and the computing time required to 
analyze determine the gas distribution within each halo 
with the full gas model (including star formation and 
feedback) within this volume. 

We proceed in the following manner. First, we mea- 
sure the mean (l&cfc) and standard deviation dbck of the 
background integrated SZE within a circular aperture of 
radius 9 centred around each cluster in a sample of im- 
ages generated from the lightcone. We do this for a range 
of angular radii, varying 9 from 0.5 to 15 arcminutes. We 
then obtain the total scatter in the integrated SZE mea- 
sured within 9 for a cluster of mass M200 at a redshift z 
by combining in quadrature the scatter due to variations 
in the background flux with that due to variations in the 
internal properties of clusters, as measured in Sec. [3] 

To measure (Fbcfe) and Sbck we randomly select a sam- 
ple of 100 clusters, each with mass M200 > 10^^h~^MQ, 
over a redshift range 0.2 < z < 1 from the lightcone. For 
each cluster we generate two images, one including and 
one omitting all the foreground and background clusters 
within a one degree radius around the cluster. We then 
remove the latter from the former, leaving just an im- 
age of the interlopers. We wish to investigate the impact 
of the SZ background due to dim clusters along the line 
of sight, therefore we also remove all clusters with an SZ 
signal to noise > 5 from the images, where Y is measured 
within 9{R2oo,z) and we assume an instrument noise of 
10/iK. These clusters are then removed in the following 
manner. First, the central pixel of the interloper is iden- 
tified and the flux ymax recorded. Next, a beta model of 
the form 

y{9)^y^aM + 9y9l)-"-^ (9) 

is deducted from the image, where 9c ~ Ryir/20Da, and 
Rvir and the angular diameter distance. Da, are taken 
directly from our lightcone halo catalog (we remove the 
cluster signal out to the angular size corresponding to 
the virial radius of each cluster) . 
Note that we do not add the primary CMB signal or 
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Fig. 4. — The mean ((Ybck): crosses) and standard deviation 
{^hck I cireles) of the background SZE flux measured within an 
angular region of radius 8 centred around each bright cluster in 
our sample of images, due to dim, unresolvable clusters (see text). 
Solid and dashed lines give power-law fits to the results, as given 
in Eqns. [TOlandlTTI 

other sources of 'noise' - radio point sources, instrumen- 
tal noise, galactic dust emission and atmo spheric emis- 
sion (see, for example, iSehgal et all 120071 ) - into these 
images. In a subsequent paper (Shaw & Holder, in prepa- 
ration) we will investigate in more depth t he accuracy 
with which matched-filtering schemes (e.g. iMelin et al.l 
can measure clusters masses, incorporating some 
of these effects. Here we focus solely on sources of scatter 
in the Y — M relation that are due to variations in the 
internal properties and spatial distributions of clusters, 
and thus cannot be suppressed. 

We measure the mean and standard deviation of the 
background integrated SZE at each value of 9 over the 
selected cluster sample. The results are displayed in Fig. 
m The crosses represent the mean integrated flux (Ybck) 
within 9. The circles represent the standard deviation 
Shck around the mean for each 9. We find that both are 
well described by a power-law, with 



and 



(Ybck) = lO-^^^^fli*^ (10) 

<5bck = 1O-5-401-61 , (11) 

plotted as solid and dashed lines, respectively. We note 
that both (Ybck) and 5bck are dependent on the matter 
power spectrum a nd thus on the value s of ^Im and as (as 
demonstrated by iHolder et al.l 12007) and therefore the 
values measured here are only relevant for the cosmology 
assumed by our simulation. 

Using Equations [10] and [11] and the fitted values for 
the Ya — -^200 scaling relations measured in Section. [3l 
we are now able to calculate the mean total flux that 
would be measured within an angular radius 9{Ay,z) 
for a cluster of mass M200 at redshift z: 

Ytot{M2oo, 9) = (ncfe(0)) + rclus(Af200, 0, z) , (12) 
where 

J^clusijW200,t/, — — -J—— iU 

aj^yz) 



1014/i-iMf, 



and = 0(A, z) is the angular radius corresponding to 
a cluster with radius i?A at angular diameter distance 
dA{z) (note the units of Fcius are arcmin^, see Eqn. [3]). 
A(A) and a(A) are the corresponding normalisation in- 
dex and slope obtained for the Ya — -^^200 relation as mea- 
sured in Section [21 The radius at a given A is calculated 
at redshift z by relating i?A to i?200 through the NFW 
profile. We u se the mass-concentration fitting formula of 
iDolag et al.l (|2004D to calculate the mean concentration 
for a cluster of mass M200 at each cluster redshift. We 
assume that the slope and normalisation of the Y — M 
relation are independent of redshift (havin g accounted 
for the bubble scahng, E{z)). iNagail (|2006f ) have found 
this to be the case for a simulated sample of clusters 
encompassing a redshift range > z > 2. 

We can calculate the total scatter in Ytot{9, M, z) by 
combining the scatter in the mean background flux with 
the intrinsic scatter in the Y — M relation (Tym(Ay); 



atoti9,M,z) 



where 



itot 

5cliis — yclusO'YM(AY) 



(14) 
(15) 



(13) 



thus converting the fractional scatter ctym ~ <^cius/Y"cius, 
to an absolute value. 

Figure [Hi demonstrates the dependence of atot (solid 
lines) on angular radius for clusters of mass lO^'*/i~-'^M0 
and 2 x W'^h'^MQ, at redshifts z = 0.5 and 1. Also 
plotted are iJdus/Ytot (dashed) and (5bck/Ytot (dotted). 
The minimum of the dashed line corresponds to angular 
radius subtending i?5oo for each cluster mass and red- 
shift. For a 10^''/i~^Mo cluster (left panels), including 
the SZE background moves the optimal value of to a 
lower angular radii than that predicted by the intrin- 
sic cluster scatter 5c\ns alone. At z = 0.5, this decrease 
is 0.4 arcminutes, and at z = 1 it is 0.25 arcminutes. 
At higher redshifts, the cluster subtends a smaller angu- 
lar region and thus the impact of the background flux 
is lessened. This is also the case as cluster mass is in- 
creased. The right panels of Figure [5| demonstrate that, 
for a 2 X lO^^/i^^M© cluster, the background flux is small 
compared to the cluster signal, and therefore i5bck does 
not change 9 significantly. Hence, for clusters of mass 
-^-^200 > 2 X lO^*/i^^M0 variations in the SZE background 
are negligible compared to the intrinsic scatter in Y due 
to variations in the internal properties of clusters. 

Finally, in Fig. [6| we plot the optimal angle 0opt(^) 
within which the integrated SZ signal provides the least 
scatter measure of M200 as a function of redshift. The 
solid, dotted, dashed and dot-dashed lines correspond to 
1,2,5 and 10 x 10^^/i^^Mq clusters respectively. Thus 
^opt(-z — 0.5) varies between 1.5 and 4 arcminutes as 
mass increases. For z > 1, all clusters have an optimal 
angular radius less than 2 arcminutes. 

5. DISCUSSION AND CONCLUSIONS 

The Sunyaev-Zel'dovich effect is one of the most 
promising means for precise determination of cosmolog- 
ical parameters because the mapping between the inte- 
grated SZ flux and cluster mass is expected to have very 
low intrinsic scatter. Thus there is a need for concerted 
investigation into the origins and character of the scatter 
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Fig. 5. — {solid lines) The total scatter atot (Eqn. I14I I in 
the Y — M relation for Y measured within angular radius 6 for 
a 10"/i-^Mo {left panels) and 2 X lO^^h-^M© (right) cluster at 
redshifts 0.5 {upper panels) and 1 {lower). The dashed lines give 
the fractional variation in the intrinsic scatter {^dus/Ytot , Eqn. 
list , the dotted lines give the fractional variation in the SZE back- 
ground (i5bck/ Ytot), within 8 for each cluster mass. 



14 . 

12- \ 

S 10 \ 



N 



CD 

4 

2 




0.5 1 1.5 

z 

Fig. 6. — The optimal angle Sopt within which SZE flux can be 
measured for a cluster of mass lO^* {solid line), 2 X IQl* {dotted), 
5 X IQi* {dashed) and lO^^ {dot-dashed) h~^MQ to give the last 
scatter (intrinsic plus projected) in the Y — M relation cluster at 
redshift z. 

that is present in this relation, and how the scatter may 
be reduced. 

In this paper, we have presented a detailed analysis of 
the intrinsic scatter in the integrated SZE-cluster mass 
(Y — M) relation, using large simulated cluster sam- 
ples generated from a semi-analytic model of intracluster 
gas applied to the lightcone output of a high resolution 
N-body simulation, and also clusters extracted directly 
from an adiabatic SPH simulation. In particular, the 
main aims of this study were to investigate the impact 
on this relation of: 

• the choice of how cluster mass is defined and the re- 
gion within which the integrated SZE is measured; 
i.e. the impact of measuring these quantities within 



different radii i?A, where A is the overdensity rel- 
ative to the critical density, 

• incorporating energy feedback due to supernovae 
and AGN outflows into the cluster gas model, 

• variations in host halo concentration and substruc- 
ture populations, 

• the error in SZ flux measurements due to confusion 
with unrcsolvable (i.e. low Y) clusters along the 
line of sight. 

To assess the impact of the definition of Y and M, 
we measure both quantities within radii corresponding 
to a range of overdensities between A = 50 and 1500. 
We measure the slope, normalisation and scatter in the 
Ya — Ma relation using all combinations of Ay and Am 
(not just Ay = Am)- The effect of energy feedback on 
the scatter is investigated by generating three realisa- 
tions of our gas model, respectively including and omit- 
ting energy feedback and a toy model in which the gas 
is isothermal and follows the dark matter density. Each 
sample contained 1267 clusters. For comparison, we also 
take a sample of 212 clusters from the ou tput of an adia- 
batic simulation (jMuanwong et alll2001h . For each sam- 
ple, we also measure the concentration and fraction of 
mass contained in substructure of the host dark matter 
halos. 

Finally, we address the issue of the SZE background 
using the full lightcone output of our N-body simulation. 
We measure the SZE background due to low mass sys- 
tems within a range of solid angles, combining the scatter 
in the background flux with the intrinsic scatter in the 
Y — M relation (within a given angular radius) to deter- 
mine the optimal angular size within which the scatter in 
the measured Y for a cluster of given mass and redshift 
is least. Our results are summarized below. 

1. Scatter in the intrinsic Y — M relation is least when 
cluster mass is defined within as large a radius as 
possible (to a maximum of i?5o in this study) and 
the integrated SZE is measured in the range i?5oo < 
Ra < ^300- For A/50 the scatter is least for Vsoo- 
For M200, the scatter is least for Isoo- In general, 
the least scatter measure of Y is always obtained 
within a smaller radii than that at which the mass 
is defined. Wc find that the best overdensity within 
which to measure Y is given by Ay = 1.9Aj\/ + 150. 

2. Inclusion of energy feedback in our gas model 
significantly increases the intrinsic scatter in the 
Y — M relation for all combinations of Am and 
Ay. We find that this is due to significantly larger 
variations in the gas mass fraction within any ra- 
dius than in our other cluster samples. 

3. Variations in host halo concentration/central po- 
tential (for clusters of the same mass) , and the re- 
sulting impact on the cluster gas distribution, pro- 
vide a reason as to why the integrated SZE provides 
a tighter proxy for cluster mass than the central 
decrement. Constraining halo concentration (i.e. 
selecting only clusters within a narrow concentra- 
tion range) changes the Y — M relation for high 
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Ay SO that the least scatter measure of Ma is ob- 
tained by measuring the integrated SZE within a 
radius R < i?i5oo, the minimum radius probed in 
this study. 

4. Substructure increases scatter uniformly for all 
combinations of Am and Ay- Removing clusters 
with large amounts of substructure reduces scatter 
in the y - M relation by up to 20%. 

5. The mean integrated SZE around a cluster due 
to low-mass foreground and background clusters 
along the line of sight scales with angular radius 
(Ybck) io-5-66g/i-85^ ^-j^g Standard deviation 
as 6bck = 10-5-401.61. 

6. For cluster of mass M200 > 2 x lO^''/i-^M0, the op- 
timal angular radius within which one should mea- 
sure the SZE so as to obtain the tightest scaling 
with cluster mass (M200) is just the angle subtend- 
ing i?5oo at the cluster redshift (the radius at which 
intrinsic scatter due to variations in internal clus- 
ter properties is least). Below this mass, the mag- 
nitude and variance of the SZE background rela- 
tive to the cluster signal reduces the optimal an- 
gular radius (e.g by 0.4 arcminutes at z = 0.5 for 
a lO^'*/i-^M0 cluster). In addition, we have pro- 
vided a chart (Fig. [6|) giving the optimal angular 
size for measuring clusters masses using the SZE 
for a range of masses and redshifts. 

One important result of this study is that we measure 
there to be at least 20% intrinsic scatter in the Y — M 
relation due variations in the internal properties of clus- 
ters alone. When mass and SZ flux are measured at 
.R200, the scatter is just above 25%. This is higher than 
the 10-15% scatter that has been quoted by previous 
authors using small samples of clusters extracted from 
highly sophisticated simulations including complex gas 
physics (e.g. [Nagai_2006) or larger samples of less-well 
resolved clust^rslda Silva et al.| [200l iMotl et al.ll2005t 
iHallman et al.l I2007D . As noted above, we have found 
that much of the scatter we measure is due to the varia- 
tions in the gas mass fraction between clusters of similar 
mass, an effect that is greatly enhanced by the inclusion 
of energy feedback into our models. However, observa- 
tional measures of the gas fraction from X-ray studies 
of bright clusters certainly seem to suggest that there is 
indeed much variation in this quantity between clusters 
(jVikhlinin ct al. 2006; Zhang et al. 2007). Nevertheless, 
we intend to perform a rigorous comparison between the 
predictions of our gas model and the results of hydrody- 
namical simulations within the near future. 

We have shown that a significant fraction of the scat- 
ter in the Y — M scaling relation is due to variations 
in the host halo properties, specifically in concentration 



Ascasibar, Y., Sevilla, R., Yepes, G., Miiller, V., & Gottlober, S. 

2006, MNRAS, 371, 193 
Avila-Reese, V., Colm, P., Gottlober, S., Firmani, C., & 

Maulbetsch, C. 2005, ApJ, 634, 51 
Bahcall, N. & Fan, X. 1998, ApJ, 504, 1 

Balogh, M. L., Babul, A., & Patton, D. R. 1999, MNRAS, 307, 463 
Barbosa, D., Bartlett, J., Blanchard, A., & Oukbir, J. 1996, A&A, 
314, 13 



and substructure populations. Our results suggest that, 
if all halos at a given mass had the same merger history, 
then the central SZE decrement might provide the best 
measure of halo mass. However, the variations in clus- 
ter structure due to mergers greatly and preferentially 
increase the scatter in Y when measured only within the 
central regions (i?i5oo)j hence it becomes necessary to 
integrate to larger radii. We have shown that there is a 
limit to how far out one should go, however, before pro- 
jection effects due to dim, unresolvable clusters begin to 
significant effect SZE flux measurements. 

We have not investigated the impact of morphology 
and environment in this study, which may also play sig- 
nificant roles. We leave this to future work. Furthermore, 
our gas model is not able to accurately account for some 
of the complex physical processes such as shock heating 
during mergers and cooling in cluster cores, which will 
certainly introduce additional scatter in scaling relations 
jMcCarthy et al. 2004; O'Hara et al. 2006; Poole et al] 
I2OO/D . 

Finally, several studies have demonstrated there to 
be significant errors in extracting the correct value of 
Y from synthetic sky maps due to CMB confusion, in- 
strument noise, the effect of galactic dust emission, dif- 
fuse gas outside of clusters and radio point sources, 
and systematic effects in th e cluster identification al- 
gorithm utihzed (iMelin et al. 2006; Hallman et al. 200'^; 
'Juin et al."2007"; 'S chafer et aLli2006ai ib; iPires et al.ii200a 
Sehgal et al. 2007). Combined with the intrinsic scatter 
in the Y — M as measured here, it is clear that there is 
still much work to be done to enable accurate estimation 
of cluster masses and thus achieve the tight cosmological 
constraints that are envisaged from SZ cluster surveys. 
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